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Abstract 



The crystal structures of potassium and cesium bistrifluoroacetates were determined 
at room temperature and at 20 K and 14 K, respectively, with the single crystal 
neutron diffraction technique. The crystals belong to the 12 /a and A2/a monoclinic 
space groups, respectively, and there is no visible phase transition. For both crystals, 
the trifiuoroacetate entities form dimers linked by very short hydrogen bonds lying 
across a centre of inversion. Any proton disorder or double minimum potential can 
be rejected. The inelastic neutron scattering spectral profiles in the OH stretching 
region between 500 and 1000 cm^^ previously published [Fillaux and Tomkinson, 
Chem. Phys. 158 (1991) 113] are reanalyzed. The best fitting potential has the 
major characteristics already reported for potassium hydrogen maleate [Fillaux et 
al. Chem. Phys. 244 (1999) 387]. It is composed of a narrow well containing the 
ground state and a shallow upper part corresponding to dissociation of the hydrogen 
bond. 
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Introduction 



In most of hydrogen bonds AH- ■ -B, it is widely accepted that the potential 
energy for the motion of the hydrogen atom has two minima [1-5]. For symmet- 
ric systems A- • -H- • -A, the two wells can be equivalent and proton tunnelling 
may occur. For the shortest hydrogen bonds the distance between the two 
minima decreases and the potential barrier between them may disappear [6]. 
There remains a single well at the centre and such symmetric structures can 
be regarded as "intermediate-states" for proton transfer in chemical reactions 
or biological processes [7,8]. Our motivation for the studies of the title com- 
pounds reported below is to shed light on the potential function experienced 
by protons involved in strong symmetric hydrogen bonds. 

A classic example of a symmetric intramolecular hydrogen bond is the mono- 
anion of maleic acid in potassium hydrogen maleate, KH(OOC-CH=CH-COO). 
This extremely short O- • -H- • -O bond, with length Rq-o — 2.437 A and with 
the proton located at the centre, has been studied with different techniques: 
X-ray diffraction [9,10]; infrared and Raman [11-15] and inelastic neutron 
scattering (INS) [16-18]. More recently, the shape of the potential function 
for the proton motion along the hydrogen bond has been determined from 
neutron diffraction and INS spectra of single crystals at low temperature [19]. 
This potential is composed of a sharp well at the centre of the hydrogen bond 
and two secondary minima at ^ 1000 cm~^ above the central minimum. These 
secondary minima are located at ^ ±0.8 A from the centre. Upon the assump- 
tion that these secondary minima occur because the hydrogen bond geometry 
is no longer stable after excitation of the u OH mode above 500 cm~^, it 
could be concluded that whilst the ground state wave function is "hydrogen- 
bonding" the excited vibrational states are "hydrogcn-nonbonding" . The po- 
tential shape is a snapshot of the proton motion during the dissociation of 
the hydrogen bond upon excitation of the OH stretching mode. The appar- 
ent paradox of a modest dissociation energy threshold observed in one of the 
"strongest" known hydrogen bonds was ascribed to partial compensation of 
the energy of formation of the hydrogen bond by a strain energy arising from 
the planarity of the maleate ring. This potential function was criticized by 
Wilson et al. [20] on the basis of plane-wave DFT calculations. However, such 
calculations performed over a very limited range of proton displacement, cor- 
responding to potential energy far below the ground state, are not conclusive 
[21]. 

So far, potassium hydrogen maleate is the only reported example we are aware 
of this unforeseen potential shape and whether this system is unique or if 
similar potentials also apply to other strong symmetric hydrogen bonds is an 
open question. 
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By contrast, a totally different potential shape has been proposed for potas- 
sium and cesium hydrogen bistrifluoroacetates, KH(CF3COO)2 and 
CsH(CF3COO)2, which are also classic examples of strong and symmetric hy- 
drogen bonds at room temperature [6] with lengths Rq-o — 2.435 and 2.38 A, 
respectively [22] . The vibrational spectra of these salts have been thoroughly 
investigated [23,24]. Upon comparison of infrared, Raman and INS spectra, 
the vibrational modes for protons were assigned at ^ 800 cm^^ {v OH), 1300 
cm~^ (7 OH) and 1600 cm~^ (5 OH). For both systems, the v OH band profiles 
observed with INS and in the infrared are almost identical. They show complex 
multi-component profiles with main components at 700, 800 and 880 cm~^ 
and weaker components in the range 500-1000 cm~^. These profiles were re- 
garded as mainly due to Evans-type [25,26] or Fano-type [27] interactions with 
internal modes of the CF3COO entities. For the Cs derivative, a very sharp 
INS band at ~ 90 cm~^ was tentatively assigned to proton tunnelling along 
the OH stretching coordinate and an asymmetric double minimum potential 
was proposed, quite at variance to the case of potassium hydrogen maleate. 

However, with hindsight, such a double well is rather unlikely for very short 
hydrogen bonds because the minima located at ^ ±0.5 A off centre (see Fig. 7 
in Ref. [24]) are quite incompatible with the Rq-o length. The O — H distance 
less than 0.7 A is quite unrealistic. In addition, there is no evidence of v 
OH tunnelling transition for the K derivative, whereas band profiles, crystal 
structures and hydrogen bond lengths are very similar. 

The localization of the proton in an asymmetric double well is also in con- 
flict with the known crystal structures determined with X-ray diffraction [22]: 
hydrogen bonds are centrosymmetric for the two salts at room temperature. 
However, crystal structures at low temperature are unknown and a symmetry 
breakdown cannot be excluded. In order to determine the proton location we 
have used the single crystal neutron diffraction technique. The crystal struc- 
tures of the two salts at 20 K or 14 K and at 298 K (room temperature) 
are presented below. In both systems, hydrogen bonds are symmetrical at 
any temperature, with protons clearly located at the centre. It is shown that 
potential functions largely inspired by that of potassium hydrogen maleate, 
though significantly different, account for the OH stretching INS band profiles. 



1 Crystal structures 

Single crystals were obtained by slow recrystallization from aqueous solutions. 

These very hygroscopic colorless crystals were handled under dry atmosphere. 
For neutron diffraction measurements, approximately cubic samples (3x3x3 
mm^) were cut from large crystals and tested at room temperature. Samples 
were loaded into aluminum containers that were then mounted in a cryostat. 
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For measurements at low temperature, they were cooled down with a flow of 
helium vapor. 



Measurements (see Table 1) were carried out on a Stoe four-circle diffractome- 
ter 5C2 at the Orphee reactor (Laboratoire Leon-Brillouin) [28] . The incident 
neutron wavelength was A = 0.8305 A. Absorption corrections were ignored. 
Data analysis was carried out with the computer package CRYSTALS [29] 



Crystal structures at any temperature are similar to those previously deter- 
mined with X-ray diffraction [22]. The potassium and cesium salts belong 
to the 12 /a and A2/a monoclinic space groups, respectively, both with four 
entities in the unit cells (see Figs. 1 and 2). The two acetate residues are 
crystallographically equivalent for both salts. They are linked by short hydro- 
gen bonds lying across a centre of inversion (see Figs. 3 and 4). Hydrogen 
bond lengths are identical at low temperature: Rq-.-o = 2.436(3) and 2.436 (4) 
A, for the K and Cs derivatives, respectively (see Tables 3 and 6). At room 
temperature, they are slightly different: 2.432(3) and 2.444(2) A, respectively. 
However, the difference is tiny compared to the thermal factors for O atoms 
(fa 0.05 A2, see Tables 4 and 7). 



At any temperature there is no evidence for splitting of the proton proba- 
bility density that could be attributed to disorder over different sites or to a 
double well potential. At low temperature, the temperature factors for oxygen 
atoms are much smaller than those for protons. As the latter compare to the 
mean square amplitudes of displacements for proton oscillators with averaged 
frequencies of 1000 cm~^ (see Tables 4 and 7), they are primarily due to 
proton dynamics in a single well. At room temperature, the temperature fac- 
tors increase substantially (see Figs. 3 and 4). For protons these factors can 
be regarded as combinations of the proton contribution, virtually unchanged 
compared to that at low temperature, with those of the oxygen atoms (those 
labelled 0(2) in Tables 4 and 7), which increase dramatically. Therefore, even 
at room temperature, the temperature factors are not compatible with any 
realistic double well potential, with minima separated by more than 0.5 A 
[24]. 



Finally, trifiuoroacctate entities have different conformations for the two salts 
(see Figs. 3 and 4). A C — F bond is trans with respect to the CO(H) bond for 
the K salt. One of the C — F bonds is virtually perpendicular to the carboxylic 
plane for the Cs derivative. 
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2 INS band profiles and proton dynamics 

The (CFaCOO- • -H- • -OOCCFa)- entities occupying Ci sites, the OH stretch- 
ing and bending modes are only infrared active (^41* + Bu) [23] . The stretching 

mode gives a broad continuum of intensity centered at ~ 800 cm"-*-, whereas 
bending modes are extremely weak and barely visible. This assignment scheme 
is corroborated by the correlation between i> OH frequencies and hydrogen 
bond lengths, Rq-o [6]- It is also consistent with the OH stretching at 600 
cm~^ for the maleate salt [19]. 

In the crystal structures, the shortest distances between protons are 4.34 and 
6.72 for K and Cs salts, respectively. These protons are further isolated by 
the stacking of CF3 entities and K"*" or Cs+ ions. Dynamical coupling between 
protons is neghgible. The great similarity of the u OH profiles observed in 
the infrared and with INS (free of symmetry related selection rule) confirms 
that crystal field splitting is not resolved [24]. Moreover, the profiles are al- 
most identical for the two salts (see below), though they have different crystal 
structures. We can safely conclude that the Au — Bu splitting is negligible 
compared to the intrinsic width of the u OH profile. 

The INS spectra in the OH stretching region, from Ref. [24], are presented in 
Fig. 5. The raw data have been reanalyzed with a better background treat- 
ment. The spectra, renormalized with respect to the amount of sample in the 
beam, show very similar intensities, as anticipated for these closely related 
salts. The OH stretching profiles were decomposed into gaussian components 
(see Tables 8 and 9) and a linear baseline was subtracted. The grey compo- 
nents correspond to Raman bands previously attributed to 6 CF3 modes [24]. 
The remaining bands, for example at ^ 602 (597), 704 (704), 799 (798), 871 
(880) and 955 (950) cm"^ for the K (Cs) derivative, are attributed to the OH 
stretching modes. 

As anticipated for hydrogen bonds with equal lengths, the u OH profiles are 
virtually identical. Only in the 820-920 cm^^ range, two components at 850 
and 870 cm~^ are clearly visible for the Cs derivative, while there is only one 
component for the K analogue. As rather strong Raman bands were observed 
at 850 and 854 cm~^ for the potassium salt [24], we conclude that the INS 
component at 870 cm^^ for the K derivative is an unresolved comprise of S 
CF3 and u OH modes (see Table 8 and 9). It sounds logical that different 
frequencies for the S CF3 modes arise from different conformations of the CF3 
groups (compare Figs. 3 and 4). 

The rather complex u OH profiles are thus supposed to arise primarily from 
the internal dynamics of the hydrogen bond. At first sight, the slowly varying 
spacing of the components, namely 102 (107), 95 (94), 72 (82), 84 (70) 
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cm"-'^ in Tables 8 and 9, suggests a Franck-Condon-like progression due to 
strong coupling with O- ■ -O modes [30-33]. However, the dynamical separa- 
tion, analogous to the Born-Oppenheimer approximation for electrons and 
nuclei, is irrelevant for strong hydrogen bonds with OH stretching at low fre- 
quency. Therefore, by analogy with the seminal case of potassium hydrogen 
maleate [19], we assume that the effective potential function for the proton 
can be decomposed into a deep and narrow well containing the "bonding" 
ground state and a shallow quasi harmonic potential for "non-bonding" ex- 
cited states. Then, the relative intensities of the components can be ascribed 
to a substantial shift of the minimum of the upper shallow curve with respect 
to the centre of the hydrogen bond. With this model, the best fitting poten- 
tial function to the observed spectral profiles of the K derivative is presented 
in Fig. 6 and Table 10. The fitting potential for the Cs derivative should be 
virtually the same. The details of the numerical calculations are given in the 
Appendix (see below). 

We have limited the number of adjustable parameters in the potential func- 
tion to the number of observed transitions. Since relative intensities were also 
included in the model, the number of experimental values exceeds the number 
of parameters. The potential function was decomposed into a narrow gaus- 
sian well with a depth of ~ 1000 cm~-^ and a rather flat harmonic potential 
whose fundamental frequency corresponds to the mean level-spacing. The lin- 
ear term displacing the upper wave functions with respect to the ground state 
allows for adjustment of relative intensities. The small cubic term accounts 
for deviations from harmonicity of the upper potential. 

The maximum difference between observed and calculated frequencies in Table 
10 is Az//z/ ^ 2%. The agreement between observed and calculated intensi- 
ties is less satisfactory, with deviations as large as ~ 30%. These differences 
emphasize the limits of the model that ignores the complex interactions be- 
tween u OH and 6 CF3 modes. Let us recall that the incoherent scattering 
cross section of fluorine atoms is negligible and the coherent cross section is 
rather modest (~ 4 barns), compared to the incoherent cross section for H 
atoms (~ 80 barns) [34]. Therefore, INS intensities for the CF3 modes are 
almost totally borrowed from the OH mode, thanks to vibrational coupling 
[24] . Consequently, the estimated frequencies and intensities obtained by sim- 
ple decomposition into gaussian components presented in Tables 8 and 9 are 
different from what they should be in the absence of such interactions. We 
did not try to model this complex interaction. Nevertheless, we consider that 
coupling with CF3 modes should not alter the main features of the potential 
function. 

The potential function represented in Fig. 6 is totally unconventional and we 
are not aware of any other reported example with similar shape. Although 
this is a particular fit with a chosen analytical function, we consider that any 
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alterative parameterization should yield basically the same overall shape for 
the best fit to the observed data. 

The potential asymmetry is not in conflict with the centrosymmetric dimers in 
the crystal structure. At low temperature, only the ground state is populated. 
The minimum is located at the center of the hydrogen bond and the wave 
function is virtually symmetrical. At room temperature, the population of 
excited states 5%) is too small to be distinguished with neutron diffraction. 
The crystal structure averaged over space and time remains centrosymmetric 
and the tiny contribution of excited states is embedded in the temperature 
factors. 

In the ground state, the proton is well localized at the centre of the hydrogen 
bond. The estimated vibrational mean square amplitude of 0.02 com- 
pares quite favorably to the temperature factors given in Tables 4 and 7. In 
the excited states, proton transfer to one of the oxygen atoms breaks down the 
symmetry of the hydrogen bond. The mean position at ^ 1 A from the centre 
is clearly incompatible with the hydrogen bond length Rq - o ~ 2.4 A when 
the proton is in the ground state. Such a large displacement of the hydrogen 
atom is possible if the hydrogen bond is broken in the excited states and the 
spatial extension of the wave functions over several A, greater than the O- • -O 
distance at equilibrium, is representative of the dissociation of the salt. 

The hydrogen bistrifiuoroacetate complex in the OH excited state is tenta- 
tively represented in Fig. 7. The proton is transferred to one of the two en- 
tities. The extension of the wave functions suggests that there is a complex 
combination of translation along the hydrogen bond and rotation around the 
C — OH bond. At the same time, there may be some twisting around the C — C 
bond to keep the proton close to the O- ■ -O direction. These complex dynam- 
ics are possible because the upper level spacing of 100 cm~^ is comparable 
to the frequency of the torsional mode [24] . Proton dynamics on the one hand, 
internal and O- ■ -H- ■ -O modes at low frequency, on the other, are on the same 
timescale and can interact strongly. 

The asymmetry of the potential in Fig. 6 contrasts to the symmetric potential 
for potassium hydrogen maleate [19]. Presumably, the electronic structure of 

the hydrogen maleate ring preserves the symmetry in the excited states. In 
contrast to the naive representation of the dissociated excited state proposed 
in Fig 10 of Ref. [19], the excited state should be better regarded as a quantum 
superposition of the two configurations corresponding to proton transferred to 
one oxygen atom or the other. 

The effective potential functions discussed above were determined upon the 
assumption of a bare proton, m = 1 amu, while coupling with heavy atoms 
were ignored. This approximation is reasonable only for the ground state. For 
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upper states, we should include dynamical coupling with internal degrees of 
freedom of the trifiuoroacetate entities, according to the scheme in Fig. 7. 
Then, the effective oscillator mass should be significantly greater than 1 amu, 
the effective upper potential should be steeper than shown in Fig. 6 and the 
spatial extension of the wave functions should be diminished. The modelling 
of the effective mass requires more information than available at the present 
stage of investigation. However, we suspect the overall shape should remain 
qualitatively the same. 

Our conclusions are quite at variance from those proposed in a recent study 
of the crystal structure of 4-cyano-2,2,6,6-tetramethyl-3,5-heptanedione, a new 
example of a very short symmetric intramolecular hydrogen bond with Ro - o = 
2.393 A at 100 K [35]. The assignment of the v OH mode to an INS band 
at 370 cm~^ is sustained by density functional calculations and magic angle 
spinning (MAS) solid-state NMR. The authors conclude (surprisingly) that 
the potential for the hydrogen bonding proton should be of the low-barrier 
double minimum type, "Although the neutron [diffraction] data better fit a 
single anisotropic thermal ellipsoid..." Unfortunately, there is no quantitative 
information on the thermal factors of the proton for comparison with other 
systems and visual examination does not reveal any significant difference com- 
pared to maleate or trifluoroacetates. Clearly, there is not enough information 
to conclude on the shape of the potential. 



Conclusion 

The crystal structures of K and Cs hydrogen bistrifiuoroacetates determined 
with the single crystal neutron diffraction technique confirm that the very 
short hydrogen bonds linking trifiuoroacetate entities are symmetrical at room 
temperature and at 20 K or 14 K. There is no evidence for proton disorder 
or phase transition. Any double minimum potential for the proton can be 
rejected [24]. 

The complex INS spectral profiles for the OH stretching modes are tenta- 
tively decomposed into gaussian components. Comparison with Raman spec- 
tra allows us to distinguish bands arising from the b CF3 modes in the same 
frequency range. We can thus assign several OH stretching transitions and 
propose a best fitting potential largely inspired by that previously determined 
for potassium hydrogen maleate [19]. This potential is composed of a narrow 
well at the hydrogen bond centre containing the bonding ground state and 
a rather flat quasiharmonic upper part corresponding to non bonding states. 
The rather low energy dissociation of this strong hydrogen bond is similar to 
that observed previously for potassium hydrogen maleate. Apparently, it is a 
property of this type of hydrogen bonding, rather than a consequence of the 
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particular structure of the maleate ring. 

In contrast to potassium hydrogen maleate, symmetry breakdown occurs for 
hydrogen bistrifluoroacetates upon excitation of the v OH modes. We conclude 
that the proton is transferred to one oxygen atom or the other. Presumably, the 

preserved potential symmetry of the hydrogen maleate in the excited states 
is a consequence of the particular structure of the ring. The excited states 
should be regarded as quantum superposition of proton transferred to the two 
oxygen atoms of the hydrogen bond. 

We conclude that the concept of hydrogen bonding for the ground state and 
hydrogen nonbonding for excited states could be of general relevance for strong 
symmetric hydrogen bonds. 



Appendix: Calculation of energy levels, wave functions and INS in- 
tensities 

The variational method is the most appropriate to determine analytical po- 
tential functions fitting any given energy level scheme. The expansion of the 
eigen functions with harmonic basis sets allows us to calculate all matrix ele- 
ments of interest for vibrational spectroscopy. In this appendix we gather the 
formulae to resolve the Schrodinger equation for a dimensionless particle with 
effective mass m* experiencing a potential V{x) along the x coordinate: 

d^^ 

with 

V{x) - V^{x) + Vg{x) ^j^vix' + j^ai exp(-6, x^) , (2) 

1=1 1=1 

where Vp is the polynomial potential and Vq the sum of three Gaussians, to 
allow for a rich variety of potential forms. We also propose a straightforward 
method to test the accuracy of the eigenf unctions. 

In order to construct a basis set, let uoq/{2tx) be a frequency and consider 
the dimensionless variables ^ = ax, with a = ^mujo/h, V = V/{huJo) and 
E = E/{hujQ). Then 

1=1 
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with 



Vi = Vi/{hlxlQQ^) 



Here, the scaled energies are such that the difference between energy levels 
of the harmonic oscillator is 1 instead of 4 in Ref. [36]. Our constants vi, 
divided by a factor of 2 compared to [36] are such that the potential energy 
of the harmonic oscillator V — ^ jl corresponds to vi — l/25i^2 (where 5i^2 is 
the Kronecker symbol). This choice of scaling sounds more natural and clear. 
Then, a = ao ^Jv^nf can be expressed in units with cto = 0.17273 and 
= ujq/{2ttc) in cm~^ units. 

In the new variables Eq. (1) reads: 

Id^^ ~ 



with %l) — ^Ja ^ normalized with respect to ^. 

For a harmonic oscillator the normalized solutions of Eq. 3 are: 

u^{i) = (2"n!V^)-^/2 En{i) exp(-eV2) n = 0, 1, . . . , (4) 



{iir„} being the Hermite polynomials. The functions form a suitable 

orthonormal basis set for localized solutions For numerical calculations the 
dimension N of the basis has to be finite. We found A?^ = 60 is suitable to 
calculate the 10 lowest energy levels with good accuracy (see below). The 
matrix elements of h in Eq. 3 are: 



hn,m = {n\h\m) 



The matrix elements used to calculate the eigenstates for a given oscillator 
mass and a sixth order polynomial potential or a Gaussian potential can be 
found in Refs. [36] and [37], respectively. As these papers may be nowadays 
difficult to obtain, we present the algorithms in condensed form, with an al- 
ternative definition of the scaling factors for the polynomial form. 
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For the polynomial potential V — Vp we obtain: 

Kn = {n + l){l + d2) + l{2n^ + 2n+l)v^ 
+ J(20n3 + 30n2 + 40n + 15){;6 ; 
= (Ci + ^nv3 + l{2n^ + 1) v^) ; 

hl,n-2 = (n - 1) ( - ^ + ^2 + (2n - 1){;4 

+ ^ (n^ - n - 1) {^e) ; (5) 

Kn-s = (^3 + i - 1) ^a) ; 

^^,n-4 = i v^n(n-l) (n-2) (n - 3) (i;4 + | (2 n - 3) v^) ; 
_ 1 / n (n - 1) {n - 2) (n - 3) {n~^ . 

'^n.n-S — 5 V 2 ^ ' 

K,n-G = H {n - 1) (n - 2) (n - 3) (n - 4) (n - 5) i;6 . 



The terms independent on Vi arise from the kinetic energy. All other matrix 
elements are zero except for the symmetric ones = ^n,n-/' ^ = 1, ■ ■ ■ , 6. 

For a Gaussian potential term W — a exp(— &a;^), the scaled potential isW — 

a exp(— 6^^), with a = a/{hujQ) and h = h/a^. Matrix elements Wn,m. in the 
subset {un\n=Q cire calculated according to Ref. [37]: construct a (A^, 2A^ — 1) 
auxiliary matrix G with an iterative procedure. The first row is calculated as: 



\2^{lVh) {m/2)\ 



for m(even)=0,. .., 2N-1 , (6) 



with 9 = h/{l + h). The second row is calculated as: 

Gi,m = V^Go,m~i + Vm + 1 Go,m+i ] for m (odd)=l, . . ,2N-2 (7) 



Each successive row of index n(2^n^A" — 1) depends on the two previous 
ones as: 



Gn,m = {^/mGn-l,m-l + Vm+ lGn-l,m+l " Vn- lGn-2,m); 

m = n,n + 2,n + 4, . . . , 2N — n — 1 

All elements not explicitly assigned are set to zero. Let us redefine G as the 

(A^, A^) square matrix corresponding to its first A^ columns. The procedure 
above has led to an upper triangular matrix. The elements of the lower triangle 
are obtained by symmetry G{n,m) = G{m,n), for aU n,m such that m < n. 
The (A^, A^) matrix corresponding to is = dG(b). 
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Therefore, the matrix h corresponding to the full potential V in Eq. 2 is: 

1=1 

If {En}n=o are eigenvalues in increasing order, and C is the {N, N) matrix 
whose column n is the normalized eigenvector corresponding to the eigen- 
functions are 

iV-l 

The eigenvalues and eigenfunctions in physical units are 

En = huJo En 

and 

JV-l 

= y/aJ2 Cn,mUmiax). 

m=0 

The eigenfunctions are analytical functions (although with numerically 
calculated coefficients) composed of a polynomial of order N — 1 multiplied by 
a Gaussian. Therefore, derivatives are easily obtained and the accuracy of the 
solutions can be checked by substitution in Eq. (1). As a rule of thumb, the 
last coefficients of each series {Cn,m}m=o (^ay, the last 10 for N = 60) have 
to be very small with respect to the largest of the others. This is generally 
achieved for the first few eigenfunctions to be compared with observation. 

The parameter ujq largely determines whether the truncated expansions of 
{"^n} ai'c good approximations. As a rule of thumb the exponential in Eq. 4 
should be small (say ~ e~^) at the estimated limits for the the particle position. 
This leads to lvq ~ 16 fi/(m* Ax^), where Ax is the width of the classically 
allowed region, or uq ^ 16 / {m* a'^Ax'^) with z/q and Ax in cm~^ and A units, 
respectively. For a single minimum this is obtained if huo is close to the 
first observed transition. In this case, basis sets limited to sizes N ^ 40 are 
convenient. For potentials composed of a narrow well and a shallow upper 
part, like those for strong symmetric hydrogen bonds under consideration in 
this paper, it is necessary to increase the size to iV = 60. The accuracy for 
the 10 lower eigenvalues is far beyond experimental errors and the accuracy 
of the eigenvectors is better than 1%. Further increment of N is unnecessary 
as numerical errors become larger the higher powers of x. 

The INS intensity for a transition |0) \n) at energy En is proportional to 
the scattering function 

S{Q,E) = \{0\exp{-iQx)\n)\^5{E - En). 
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For a spectrometer like TFXA [38], energy and momentum transfer, E and 
Q, respectively, are correlated as 



16.759m*' 

with E and Q in cm"-*^ and A"-*^ units, respectively. 

Therefore, for a given set of parameters {vi}, {a/}, {&/} and Uo, we are able 
to obtain the frequencies and intensities of the first few transitions with good 
accuracy. Using standard numerical methods it is possible to determine the 
values of the parameters that best fit the observed frequencies and intensities. 
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Table 1 

Neutron single crystal diffraction data and structure refinement for potassium and 
cesium hydrogen bistrifluoroacetates. A = 0.8305 A. Space groups monoclinic 12/ a 
for potassium and A2/a for cesium. Both with Z = 4. The criterion for observed 
reflections was I > 3(t(I). The variance for the last digit is given in parentheses. 



KH(CF3C00) 

20K 

a (A) 8.68(1) 
b (A) 10.023(9) 
c (A) 9.146(9) 
(3 {') 100.36(8) 
V (A3) 782.6 
(Mg m-3) 2.259 



2 CsH(CF3COO)2 

298 K 14 K 298 K 

8.78(1) 13.44(1) 13.623(8) 

10.18(1) 4.942(9) 5.033(3) 

9.28(1) 14.35(1) 14.741(6) 

99.96(9) 112.88(9) 112.46(9) 

817.0 878.4 934.0 

2.164 2.722 2.559 



Measured reflections 1979 

Independent reflections 1769 

and observed reflections 1503 

Rint 0.037 



2164 3476 3897 

1472 1977 2096 

1142 1557 1062 

0.038 0.048 0.064 



Refinement on 
R-factor 

Weighted R-factor 
Goodness of fit 
Number of reflections 
Number of parameters 

used in refinement 
Extinction coefficient 



F F 

0.040 0.045 

0.042 0.029 

1.070 1.088 

1503 1142 

76 76 

22.7(7) 17.2(8) 



F F 

0.0399 0.043 

0.0345 0.040 

1.049 1.079 

1557 1062 

76 76 

10.0(3) 11.2(9) 
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Table 2 

Atomic positions and isotropic temperature factors for KH(CF3COO)2 at 20 K 
(first lines) and 298 K (second lines). The variance for the last digit is given in 
parentheses. 





X 1 a 




z 1 c 






u.zouu 




l.UUUU 


n nn99 

U.UUZZ 




n 9^nn 
u.zouu 


n ZLRQI (0\ 


1 nnnn 

l.UUUU 


U.UoUD 






u.Doyyzi^o j 


n fii 071 ( '\\ 

U.DIU / 1 ^0 j 


u.uuzu 




n 0771 (^ \ 

U.U 1 1 1 11 J 


u.uuzool y J 


n fii 9ni (Q\ 


u.uouu 


CO) 


-U.U^-LOOl U 1 


u.uuuoy 1 (J J 


n 71 '^n^^'fi'i 

U. i lOUOl U J 


u.uuzu 




-u.U4^:oOD^^y j 




n 71 t^i 9^'fi'\ 

U. ( iOiZl^o j 


U.UZDU 


U(ij 


-u.uzooy(^ / j 


U.0oDU4(^D j 




U.UU4U 




-0.02706(13) 


0.58871(11) 


0.81966(11) 


0.0333 


0(2) 


-0.15042(7) 


0.74496(6) 


0.67180(7) 


0.0045 




-0.15011(14) 


0.74556(14) 


0.67428(14) 


0.0426 


F(l) 


0.15118(7) 


0.78039(7) 


0.61351(7) 


0.0050 




0.14415(18) 


0.78021(15) 


0.61321(18) 


0.0519 


F(2) 


0.19203(7) 


0.56895(7) 


0.65128(7) 


0.0046 




0.18745(15) 


0.57426(17) 


0.65321(17) 


0.0490 


F(3) 


0.01440(8) 


0.63619(7) 


0.47010(7) 


0.0046 




0.01325(19) 


0.63720(16) 


0.47512(13) 


0.0474 


H(l) 


-0.2500 


0.7500 


0.7500 


0.0166 




-0.2500 


0.7500 


0.7500 


0.0574 
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Table 3 

Interatomic distances in A units and angles in degrees in KH(CF3COO)2 at 20 
K (first lines) and 298 K (second lines). The variance for the last digit is given in 
parentheses. 



C(l)-C(2) 


1.5504(17) 


C(l)-F(l) 


1.3503(15) 




1.5455(19) 




1.335(2) 


C(l)-F(2) 


1.3296(16) 


C(l)-F(3) 


1.3332(17) 




1.327(2) 




1.322(2) 


C(2) 0(1) 


1.2230(15) 


C(2)-0(2) 


1.2779(16) 




1.2148(18) 




1.270(2) 


0(2)-H(l) 


1.2179(13) 


C(2)-C(l) - F(l) 


109.9(1) 




1.2159(17) 




110.08(14) 


C(2)-C(l) - F(2) 


111.9(1) 


C(2)-C(l) -F(3) 


111.57(11) 




111.51(14) 




111.87(13) 


F(l)-C(l) - F(2) 


107.86(12) 


F(l)-C(l)-F(3) 


107.08(11) 




107.78(17) 




107.18(16) 


F(2)^C(1) - F(3) 


108.29(11) 


C(l)-C(2)^0(l) 


119.4(1) 




108.24(16) 




119.61(12) 


C(l)-C(2)-0(2) 


111.4(1) 


0(l)-C(2)-0(2) 


129.21(9) 




111.36(13) 




129.03(12) 


C(2)-0(2)-H(l) 


114.2(1) 


0(2)-H(l)-0(2) 


180.0(1) 




114.87(13) 




180.0(1) 



19 



Table 4 

Thermal parameters in units for KH(CF3COO)2 at 20 K (first lines) and 298 
K (second lines). The variance for the last digit is given in parentheses. 



Atom 


Uii 


U22 


U33 


U23 


Ui3 


U12 


K(l) 


0.0027(4) 


0.0018(4) 


0.0020(4) 


0.0000 


0.0003(3) 


0.0000 




0.035(1) 


0.0263(8) 


0.0325(9) 


0.0000 


0.0122(7) 


0.0000 


C(l) 


0.00117(19) 


0.00268(19) 


0.0023(2) 


0.00024(13) 


0.00039(14) 


0.00003(13) 




0.0297(4) 


0.0329(4) 


0.0299(3) 


0.0029(3) 


0.0078(3) 


-0.0022(3) 


C(2) 


0.00157(19) 


0.00238(19) 


0.0023(2) 


0.00055(14) 


0.00060(14) 


0.00064(13) 




0.0258(3) 


0.0272(3) 


0.0249(3) 


0.0045(3) 


0.0043(2) 


0.0005(3) 


0(1) 


0.0038(2) 


0.0048(2) 


0.0035(2) 


0.00241(17) 


0.00120(17) 


0.00108(15) 




0.0360(5) 


0.0341(4) 


0.0307(4) 


0.0115(4) 


0.0088(3) 


0.0042(4) 


0(2) 


0.0033(2) 


0.0060(2) 


0.0047(2) 


0.00268(17) 


0.00194(17) 


0.00312(16) 




0.0376(5) 


0.0507(7) 


0.0418(5) 


0.0194(5) 


0.0129(4) 


0.0180(4) 


F(l) 


0.0049(2) 


0.0048(2) 


0.0054(2) 


0.00042(18) 


0.00072(18) 


-0.00222(17) 




0.0523(8) 


0.0457(7) 


0.0597(8) 


0.0074(6) 


0.0150(6) 


-0.0195(6) 


F(2) 


0.0027(2) 


0.0054(2) 


0.0057(2) 


0.00058(18) 


0.00076(18) 


0.00219(17) 




0.0353(6) 


0.0579(8) 


0.0557(7) 


0.0057(6) 


0.0135(5) 


0.0141(5) 


F(3) 


0.0054(2) 


0.0057(2) 


0.0025(2) 


-0.00028(17) 


0.00020(17) 


-0.00006(17) 




0.0608(8) 


0.0535(7) 


0.0285(5) 


-0.0032(5) 


0.0094(5) 


0.0005(6) 


H(l) 


0.0177(7) 


0.0158(8) 


0.0160(7) 


0.0034(6) 


0.0023(6) 


0.0029(6) 




0.0539(17) 


0.0577(18) 


0.0582(17) 


0.0199(15) 


0.0034(13) 


0.0152(14) 
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Table 5 

Atomic positions and isotropic temperature factors for CsH(CF3COO)2 at 14 K 
(first lines) and 298 K (second lines). The variance for the last digit is given in 
parentheses. 



iixom 


X / CL 




zjc 






-u.zouu 


n 91 1^41 (^ 

U.Z104:ll lo 1 


u.uuuu 


n nni ^ 

U.UUiO 




-u.zouu 


91 7n('zL'l 


n nnnn 
u.uuuu 


u.uo 1 / 




-U.Uloool^O ) 


n fii QT^n ^"i 


u. lozooi^o j 


u.uuoy 




-U.U±0/Ul-L^J 




n 1 i^'^Rfi('i4'\ 

U.lOOOUl 14 1 








n ^Q7^8n 111 








nil 

u.il4oo^^io ) 


U.ooUD(^0 ) 


n Ofii 1 nn 01 


U.U(D1 


b[6) 


n 1 /if?s2/i 

U. 14004(^0 j 


n 7071 sjn /11 

U. ( / ( io(^14 j 


u.ioooy^o j 


U.UUDo 




0.14651(17) 


0.7048(4) 


0.18308(15) 


0.0722 


0(1) 


0.16871(5) 


0.28677(13) 


0.08444(5) 


0.0050 




0.16517(9) 


0.2837(3) 


0.0796(1) 


0.0445 


0(2) 


-0.00501(5) 


0.18925(14) 


0.05235(5) 


0.0045 




-0.00336(8) 


0.1803(3) 


0.0543(1) 


0.0442 


C(l) 


0.07995(4) 


0.52050(11) 


0.17554(4) 


0.0031 




0.08238(8) 


0.5005(2) 


0.17453(7) 


0.0401 


C(2) 


0.08431(4) 


0.31424(11) 


0.09597(4) 


0.0027 




0.08415(6) 


0.30433(18) 


0.09498(7) 


0.0312 


H(l) 


0.0000 


0.0000 


0.0000 


0.0166 




0.0000 


0.0000 


0.0000 


0.0592 
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Table 6 

Interatomic distances in A units and angles in degrees in CsH(CF3COO)2 at 20 
K (first lines) and 298 K (second lines). The variance for the last digit is given in 
parentheses. 



C(l)-C(2) 


1.549(3) 


C(l)-F(l) 


1.331(3) 




1.5403(14) 




1.3236(19) 


C(l)-F(2) 


1.340(3) 


C(l)-F(3) 


1.331(2) 




1.326(2) 




1.324(2) 


C(2)-0(l) 


1.216(3) 


C(2)-0(2) 


1.277(3) 




1.2132(13) 




1.2746(13) 


0(2)-H(l) 


1.218(2) 


C(2)-C(l)-F(l) 


112.3(2) 




1.2221(12) 




112.3(1) 


C(2)-C(l)-F(2) 


109.20(17) 


C(2)-C(l)-F(3) 


111.56(15) 




110.24(13) 




111.63(11) 


F(l)-C(l)-F(2) 


107.6(2) 


F(l)-C(l)-F(3) 


108.10(12) 




107.61(15) 




107.41(16) 


F(2)-C(1)~F(3) 


107.9(2) 


C(l)-C(2)-^0(1) 


118.59(17) 




107.44(16) 




118.5(1) 


C(l)-C(2)-0(2) 


112.5(2) 


0(l)-C(2)-0(2) 


128.92(15) 




112.86(8) 




128.62(11) 


C(2)-0(2)-H(l) 


114.03(19) 


0(2)-H(l)-0(2) 


180.0(1) 




114.18(9) 




180.0(1) 
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Table 7 

Thermal parameters in units for CsH(CF3COO)2 at 14 K (first lines) and 298 
K (second lines). The variance for the last digit is given in parentheses. 



Atom 


Uii 


U22 


U33 


U23 


Ui3 


U12 


Cs(l) 


0.0009(3) 


0.0015(3) 


0.0019(3) 


0.0000 


0.0009(2) 


0.0000 




0.0319(6) 


0.0359(7) 


0.0526(9) 


0.0000 


0.0244(6) 


0.0000 


F(l) 


0.0048(2) 


0.0069(3) 


0.0061(2) 


-0.00185(18) 


0.00239(17) 


0.00167(19) 




0.0614(8) 


0.069(1) 


0.068(1) 


-0.0250(8) 


0.0275(7) 


0.0110(8) 


F(2) 


0.0095(2) 


0.0082(3) 


0.0030(2) 


0.00077(19) 


0.00246(17) 


0.0006(2) 




0.0959(13) 


0.0907(14) 


0.0369(7) 


0.0065(8) 


0.0199(8) 


0.0014(12) 


F(3) 


0.0074(2) 


0.0047(3) 


0.0068(3) 


-0.00216(19) 


0.00288(18) 


-0.00294(18) 




0.0837(12) 


0.0581(9) 


0.0777(12) 


-0.0291(9) 


0.0342(9) 


-0.0324(9) 


0(1) 


0.0038(2) 


0.0055(2) 


0.0071(2) 


-0.00185(18) 


0.00373(16) 


-0.00097(18) 




0.0344(5) 


0.0458(6) 


0.0592(7) 


-0.0061(6) 


0.0245(5) 


-0.0035(5) 


0(2) 


0.0029(2) 


0.0039(2) 


0.0068(2) 


-0.00313(18) 


0.00219(16) 


-0.00094(17) 




0.0305(4) 


0.0407(6) 


0.0634(7) 


-0.0202(6) 


0.0202(5) 


-0.0076(4) 


C(l) 


0.00369(17) 


0.0030(2) 


0.00295(18) 


-0.00084(14) 


0.00151(13) 


0.00004(15) 




0.0448(5) 


0.0407(5) 


0.0346(4) 


-0.0067(4) 


0.0148(3) 


-0.0051(4) 


C(2) 


0.00267(18) 


0.00290(19) 


0.00315(18) 


-0.00089(15) 


0.00183(13) 


-0.00025(14) 




0.0285(3) 


0.0284(3) 


0.0370(4) 


-0.0012(3) 


0.0130(3) 


-0.0011(3) 


H(l) 


0.0131(6) 


0.0176(8) 


0.0185(7) 


0.0006(7) 


0.0053(5) 


-0.0011(6) 




0.0403(13) 


0.0498(16) 


0.088(2) 


-0.0059(17) 


0.0247(14) 


-0.0011(12) 
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Table 8 

Gaussian decomposition of the INS spectrum of KH(CF3COO)2 in the OH stretch- 
ing region (see Fig. 5). 



Peak 


Width 


Gravity center 


Area 


Assignment 




cm~^ 


cm~^ 


% 




1 


12 


527 


1.9 


S CF3 


2 


36 


560 


3.4 


6 CF3 


3 


18 


602 


4.4 


OH stretching 


4 


36 


648 


5.1 


SCF3 


5 


41 


704 


17.9 


OH stretching 


6 


22 


750 


5.0 


6 CF3 


7 


20 


774 


5.0 


S CF3 


8 


38 


799 


18.1 


OH stretching 


9 


78 


871 


27.6 


OH stretching, 6 CF3 


10 


73 


955 


10.9 


OH stretching 


11 


18 


1033 


0.8 


OH stretching 



Table 9 

Gaussian decomposition of the INS spectrum of CsH(CF3COO)2 in the OH stretch- 
ing region (see Fig. 5). 



Peak 


Width 


Gravity center 


Area 


Assignment 




cm~^ 


cm~^ 


% 




1 


12 


523 


1.8 


S CF3 


2 


19 


550 


1.5 


S CF3 


3 


20 


597 


3.9 


OH stretching 


4 


26 


647 


2.6 


5 CF3 


5 


38 


704 


14.2 


OH stretching 


6 


27 


751 


4.7 


S CF3 


7 


57 


798 


30.6 


OH stretching 


8 


23 


849 


5.1 


6 CF3 


9 


53 


880 


20.1 


OH stretching 


10 


68 


950 


13.4 


OH stretching 


11 


55 


1043 


2.0 


OH stretching 
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Table 10 

Observed and calculated INS OH stretching frequencies and intensities for 
KH(CF3COO)2. The calculated ones correspond to the potential V = -185.074 a; + 
122.598x2 - 10.506x3 - 1232.04 exp(-28. 149x2), V and x in cm-^and A, respec- 
tively. For the TFXA spectrometer used in Ref . [24] , energy and momentum transfer 
are correlated as z> = Q"^ /{16.759m) with D, Q and m in cm~^, and amu units, 
respectively, an: Mrl)itrMry nulls. 





Observations 


Calculations 


Transitions 


freq. 


Q 


int. 


freq. 


Q 


int. 




(cm^^) 




(au) 


(cm^"*") 




(au) 


^ 1 


601 


6.0 


0.25 


614 


6.1 


0.38 


0^2 


703 


6.5 


0.96 


701 


6.5 


0.93 


0^3 


798 


6.9 


1.00 


787 


6.9 


1.00 


0^4 


870 


7.2 


0.76 


872 


7.2 


0.59 


0^5 


956 


7.6 


0.58 


953 


7.5 


0.37 
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Fig. 1. Projection of the crystal structure of potassium hydrogen bistrifluoroacetate 
at 20 K onto the (a, c) plane and thermal ellipsoids for atoms. 
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Fig. 3. Comparison of the hydrogen bistrifluoroacetate entities and thermal ellip- 
soids for atoms of the potassium salt at 20 and 300 K. Right: projection onto the 
mean-plane of the carboxylic entities. Left: view along the hydrogen bond direction. 
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14 K 




Fig. 4. Comparison of the hydrogen bistrifluoroacetate entities and thermal elhpsoids 
for atoms of the cesium salt at 14 and 300 K. Right: projection onto the mean-plane 
of the carboxylic entities. Left: view along the hydrogen bond direction. 
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Fig. 5. Inelastic neutron scattering spectra and band decomposition into gaussian 
profiles in the OH stretching region of KH(CF3COO)2 and CsH(CF3COO)2 at 20 
K. The filled components are attributed to other modes (see text). The residual of 
the fit is compared to the error bars. 
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-10 12 3 -1 1 2 3 



X X 

Fig. 6. Left: Calculated potential and energy levels for KH(CF3COO)2. 
V = -185.074a; + 122.598x2 - 10.506x3 - 1232.04 exp(-28.149x2), with V and 
X in cm~^ and A units, respectively. Right: The corresponding wave functions, x in 
A units. 




Fig. 7. Artistic view of the dissociation of the hydrogen bistrifluoroacetate ion in 
the OH stretching excited state 
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